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We examine the sensitiveness of the free-energy landscape of a directed polymer in random media 
with respect to various kinds of infinitesimally weak perturbation including the intriguing c ase of 



temperature- chaos. To this end, we combine the replica Bethe ansatz approach outlined in cond- 



mat/0112384, the mapping to a modified Sinai model and numerically exact calculations by the 
transfer-matrix method. Our results imply that for all the perturbations under study there is a 
slow crossover from a weakly perturbed regime where rare events take place to a strongly perturbed 
regime at larger length scales beyond the so called overlap length where typical events take place 
leading to chaos, i.e. a complete reshuffling of the free-energy landscape. Within the replica space, 
the evidence for chaos is found in the factorization of the replicated partition function induced by 
infinitesimal perturbations. This is the reflex of explicit replica symmetry breaking. 
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I. INTRODUCTION 



A very interesting problem of glassy systems with disorder and frustration is the possible instability of the glassy 
frozen states against infinitesimally weak perturbations such as an infinitesimal change of temperatures and realizations 
of quenched randomness. Such a perturbation does not bring the system out of the frozen phase but possibly changes 
the lugged landscape of the free-energy in a dramatic way. Let us call this intriguing property as fragility of the 
free-energy landscape. A class of phenomenological scaling theories started first in the context of spin-glass by Bray 
and Moore, and, Fisher and Huse ^, ^ generically implies that equilibrium states of systems with disorder and 
frustration resist against such infinitesimally weak perturbations of strength (5 ^ 1 up to a finite crossover length scale 
Lc{S) called overlap length but change into completely different states at larger length scales, resulting in the vanishing 
of the correlations between the two states. The overlap length Lc{d) diverges as 5 — > but remains finite for any non- 
zero 6. Such an anomalous response is called as chaos referring to the feature that the distance between the perturbed 
and unperturbed systems becomes infinitely large in phase space even by infinitesimally weak perturbation as the 
system size L becomes macroscopically large L/Lc{S) — > oo Q]. Unfortunately, the validity of the prediction has not 
been proven explicitly by theoretical studies except for some Migdal Kadanoff type real space renormalization-group 
(MKRG) studies |, |. Especially, the issue of temperature- chaos, i. e. the sensitivity of glassy phases with respect 
to a small change of temperature, has been of great interest because of its potential relevance for the rejuvenation 
(chaos) effects found in temperature-shift and temperature-cycling experiments ^, ^. 

The majority of the previous theoretical and numerical studies concerning the problem of the fragility of glassy 
phases has been done on Edwards-Anderson (EA) spin-glass models, which has been considered as prototypical 
model for glassy systems. While a rich amount of numerical evidences for the anomalous response to non-thermal 
perturbations has been accumulated |p, |ll|, |l2|, the intriguing problem of temperature-chaos remains very 
controversial. For the Sherrington-Kirkpatrick (SK) mean-field spin-glass model, which is the EA model embedded 
in infinite dimensional space. It is realized that saddle point solutions both with and without temperature-chaos 
p^ , p^ , pO| , pl| exist and apparently new theoretical ideas are needed. On the other hand, numerical studies report 
conflicting results [|, 0, |l|,|l3[ ||, ||, ||. 

Recently we developed an analytical scheme to study the fragility of the free-energy landscape of randomly frustrated 
systems against various kind of perturbations |Q. Especially we proposed to prove the onset of chaos in terms of 
statistical decoupling of a set of replicated partition functions, and we applied the method to the directed polymer 
in random media (DPRM). The DPRM |27) is a simple model compared to spin-glass models. In spite of this, it 
is believed to possess many of the subtle properties of glassy systems, thus, it deserves to be called as "baby spin- 
glass" [|8|. Indeed, the anomalous response of DPRM towards various kinds of weak perturbations has already been 
reported by many numerical studies ^ |3^ including a signature of temperature-chaos . DPRM belongs 

to the wide class of elastic manifolds in random media |3^, ^ |36l |37|l which encompasses a variety of physical 
systems of much interest such as the domain walls of ferromagnets with weak bond randomness and the flux 



lines in type-II superconductors with randomly distributed point like pinning centers 1 40 , 4l| , CDW and vortex lattice 
systems with weak random-periodic pinnings Q . 

The scope of this paper is to present a unified study on the fragility of the free-energy landscape of DPRM with 
respect to various weak perturbations using the replica Bethe ansatz approach outlined in , mapping to a modified 
Sinai model and numerically exact transfer matrix calculations. Our main results are the following. We find that 
infinitesimally weak perturbations amount to replica symmetry breaking terms in the effective action which lead to 
the statistical decoupling of two sets of replicas. The outcome can be naturally understood as a manifestation of 
spontaneous replica symmetry breaking following the definition of Parisi and Virasoro p7| . Interestingly enough, the 
replica approach turns out to give results quite consistent with the phenomenological scaling approach j^, |2^, [2^, |3^, 
|3l| , ^ and predict the same overlap length Lc{d). Within the replica approach, apparently different perturbations 
can be naturally classified into a few universality classes. Concerning the well known correspondence between the 
effective free-energy landscape of DPRM and the Sinai model, the statistical decoupling of replicas (chaos) naturally 
suggests the emergence of statistically independent Sinai valleys for different subsets of replicas. To examine the 
anticipated universal aspects of the anomalous response, we present and discuss the outcome of a detailed numerical 
analysis using transfer matrix methods. 

The plan of the paper is the following. In the next sections we propose a genera l fr amework to define and study 
the fragility of the free-energy landscape of randomly frustrated systems. In sectionElJ we define the DPRM model. 
In section IV, we review and summarize the previous scalin g arg uments. In section |V| we present details of the replica 
Bethe ansatz approach outlined in |26[. Then, in section VII, we present the outcom e of a n exhaustive numerical 



analysis using the transfer matrix method. Finally, we summarize our results in section VIII 
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II. STATISTICAL DECOUPLING OF REAL REPLICAS 



In this section we discuss a general strategy to study the sensitivity of the free-energy landscape of a generic class 
of systems. The free-energy _F of a random system is a random quantity with certain mean and variance. Let us 
denote the deviation of the free-energy of a given sample from the mean as, 

AF = F-F. (1) 



Here and hereafter (• • •) denotes the average over different realizations of randomness. 

Now let us consider two systems say A and B. Initially they are prepared as two identical copies with the same 
randomness, temperature and other parameters. Such systems are called as real replicas. We are interested in how 
the statistical correlation between A and B changes by introducing a perturbation of strength S. Then it is useful to 
define a disorder-averaged correlation function. 



/AFl(i)^A^^|(i) 
If the correlation function vanishes at large length scales 

lim lim Ci.(L,(5) (3) 

it implies the free-energy landscape of A and B decorrelates completely. If the statistical decoupling between A and 
B happens even with arbitrarily weak perturbation (5 <C 1, we say that there is chaos. 

Now let us consider an equivalent definition of chaos which is more suited for analytical approaches based on the 
replica method. Let us suppose that each of the systems A and B are replicated further into n replicas and consider 
the disorder-averaged partition function of the total system Z'^j^^{L). As noticed by Kardar [Q, if an analytical 
continuation for n ^ is possible, the disorder average of such a partition function can be identified as the generator 
of cumulant correlation functions of sample-to-sample fluctuations of free-energies |Q . Thus the complete knowledge 
of the disorder average of the replicated partition function allows one to obtain the distribution function of sample-to- 
sample fluctuation of the free-energy p5| . In our present context, the disorder-averaged partition function Z\^g(L) 
generates cumulant correlation functions of the total free-energy as the following. 



limlnZ^+5(L) = n\uZA+B{L) + —[\nZA+B{L)]l + ... + ^[\nZA+B{L)fc 

^2 



n- 



P. 



+ ... + —[^l3AFA-pBFBfc + ... (4) 

pi 

where [. . stands for p-th cumulant correlation functions of the total free-energies —PaFa — PbFb with Fa{L) and 
Fb{L) being free-energies of subsystems A and B respectively and (3a and (3b being inverse temperatures of A and B 
respectively. 

Obviously, the decorrelation of the free-energy fluctuations between A and B is equivalent to the factorization of 
the replicated partition function. 



hm hm hm Z^+^(i,<5)- x Zl{L) (5) 

Note that if the latter result holds, automatically(||) holds too. An important remark is that the order in which limits 
are taken is crucial to obtain sensible results: the limit n must be taken before than the thermodynamic limit 
L — > oo and finally the limit ^ — > must be taken. In what follows, we will use dq) as our definition of chaos in the 
replica approach. We have to stress, though, that this definition is general and holds for generic random systems. 

The above definition of chaos implies that it can be regarded as a spontaneous symmetry breaking phenomenon. 
If the perturbation is absent, A and B are equivalent and one expects the exchange symmetry A <-> _B to be present. 
One also expects to have permutation symmetry among the replicas associated to each group A or B. Such an 
invariance under permutations is usually called in short as replica symmetry. However, in general, it turns out that 
the disorder-averaged replicated partition function of the 2 x n replicas Z^^g without any perturbation has an even 
higher symmetry: it is invariant under any permutation among the 2 x n replicas. Now, if (||) holds, this higher 
symmetry is reduced: after having introduced a perturbation the permutation symmetry remains at most within each 
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subset associated with A and B. Thus in order the this phenomena happens, the perturbation should show up in the 
repHcated partition function as a symmetry breaking term which tries to break the fuU permutation symmetry. Now, 
the definition (^) tells us that this symmetry breaking happens even with an arbitrary weak perturbation. Therefore 
chaos defined as (|^) is a spontaneous replica symmetry breaking phenomenon. We note that such a definition of 
replica symmetry breaking was introduced under the name explicit replica symmetry first by Parisi and Virasoro [4^ 
who tried to give a sound thermodynamic definition for the replica symmetry breaking phenomena known in the 
saddle point solutions of mean-field models M, of a class of glassy systems. 



III. MODEL 



We study DPRM in 1 + 1 dimensions which is described by the following Hamiltonian in the continuous limit, 

Ho[V,h,^] ^ [ dz 
Jo 

The scalar field (j) represents the displacement of the elastic object at point z in a f -dimensional internal space of size 
L. We assume that the field is a single-valued function of z which means that oriented objects with no overhangs 
are considered. In the following, we assume that one end of the string is fixed as 0(0) = while the other end (j>{L) 
is allowed to move freely. The 1st term in the Hamiltonian is the elastic energy, k being the elastic constant. The 
random pinning media is modeled by the quenched random potential Vo((/), z) with zero mean and short-ranged spatial 
correlation, 

Vo{^ = Voicj), z)Vo{cb\ z') = 2D5{,p - (j,')5{z - z') (7) 

Many exact properties of this 1 + 1 dimensional model are known [ p7) . It is in the frozen phase at all finite temperatures 
in the sense that its scaling properties are always governed by the T = glassy fixed point. 

We implement the basic strategy explained in the previous section as the following. First, we start with a system 
of two real replicas, say A and B, whose configurations 4>a{z) and 4>b{z) are subjected to exactly the same random 
potential and temperature. Second, we apply small perturbations to them. In the present paper we consider 5 different 
kinds of perturbations: 

I) Tilt field ll^: A and B replicas are subjected to a tilting field of opposite sign —h(l)A[L) -f h(f>B{L) with ft, <C 1. 

II) Explicit short-ranged repulsive coupling | |2^ , Q: A and B replicas are subjected to explicit repulsive short-ranged 
interaction e J^^ dzS((j)A{z) — 0b(z)) with < e ^ 1. 

Ill-i) Decorrelation of random potential p9| , the random potential of B is made from that of A as Vb — {Va + 
5V')/w\ + 5^ where \5\ <C 1 and V' follows the same Gaussian distribution as V. Then Vg(0, z) = and 
Vg(/',z)Vg'(0',z') = 2DGG'5{(ly - (t)')5{z - z') with Daa = Dbb = D and Dab = D/VTTP < D. 

Ill-ii) Random Tilt Field: A and B are subjected to statistically independent weak random tilt field. 

Ill-iii) Temperature difference |^ : slightly different temperatures Ta = T + ST Tb = T — ST for A and B respectively 
with ST/T < 1 



d(j){z] 
dz 



VoWz),z) 



(6) 



IV. DROPLET SCALING APPROACH 

We first review and discuss the scaling approach picture j|, |9[ |2| for the problem of the anomalous 

response. Let us consider a simple-minded picture consisting in the deepest valley corresponding to the ground state 
configuration and many branched valleys of low (free-)energy excitations which for given longitudinal size L differ 
from the 'ground state' over a transverse size uo{L/Lq)'^, C being the so called roughness exponent. Note that we 
have introduced a characteristic longitudinal length Lq, which should be understood as the Larkin length |]50| beyond 
which pinning becomes important, as well as its associated transverse length scale uq- The free-energy gap of these 
excited states with respect to the 'ground state' is expected to scale typically as, 

Af'f P = UoiL/Lof. (8) 

Here Uq is the energy scale associated with the Larkin length and 6 is the stiffness exponent which is related to the 
roughness exponent C by the exact scaling relation. 
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e = 2C-i. (9) 

In a 1 + 1 dimensional system these exponents are believed to be exactly = 1/3 and C = 2/3 ^ The 
probability distribution function of the free-energy gap A_Fi is expected to have a natural scaling form, 

..,AnMAf.,=,(j^)ji|^, ,10, 

with non-vanishing amplitude at the origin, 

m > 0, (11) 

which allows rare, gap- less excited states |^ psf . 

Let us now consider a generic perturbation which triggers an excitation from the 'ground state' with a free- energy 
gain of order, 

su(^y, (12) 



L 



in the infinitesimally weak perturbation limits 



5U/Uq 0. (13) 



In the following we consider perturbations such that a > 9. Under the influence of such a perturbation, the system 
in the deepest valley may jump into other valleys with free-energy gap A^F if the possible gain of free-energy due to 
perturbation ( p^ becomes larger than the original free-energy gap itself. The probability of such an event is estimated 
as, 



5U(L/Lo) 



Pjump(i,<^C/)= / pL{^FL)d^F^[j-j^^j , (14) 







L V 



with a characteristic length scale called overlap length, 



L,{dU) ^ Lo i^— j as SU/Uo^O. (15) 

Let us also define a characteristic transverse length scale which is conjugate to Lc{SU), 

«c(<5C/)=-o(^^)'. (16) 

It is important to note that the above expressions make sense only for short enough length scales L ^ Lc{SU). 
In this regime the effect of the jumps on physical quantities can be analyzed in a perturbative way because the 
probability of a jump is small enough. Let us call this regime weakly perturbed regime. However, in strongly perturbed 
regime L 3> Lc{dU), perturbative treatments will fail because jump events will happen with probability one. The latter 
implies that after having applied the perturbation the free-energy landscape is drastically different from the original on 
lengthscales larger than the overlap length. The overlap length ( p^ ) diverges as SU /Uq — > with exponent —l/{a—9) 
which is sometimes called as chaos exponent []67t but remains finite for arbitrary small strength of perturbation 6. 

A. Uniform Tilt Field 

We first consider the application of a uniform tilt field h to the end-point of the real replica B at z = L by which 
the statistical rotational symmetry is violated. In the presence of the tilt field the Hamiltonian becomes, 

Ha+b = Ho[Vo, c^a] + Ho[Vo, M - Km [ dz^^§^ (17) 

. n dz 
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The unperturbed Hamiltonian Hq is given in (|6|). If the string makes a jump responding to the uniform tilt field over 
a distance of order uq{L / Lq)'^ into the next valley, it obtains an energy gain of order /iu„iMo(L/-Lo)^- Thus the unit 
for the gain in energy (|l2|) reads as 6U = h^niUo with characteristic exponent a — — 2/3. Therefore we find the 
overlap length ( |l5| ) to be, 

ic(/iu„i) ~ io (^^) ' (18) 

This result was previously obtained by Mezard in by using essentially the same argument and supporting his 
result by a numerical transfer matrix calculation. 



B. Explicit Repulsive Coupling 

The other perturbations that we consider do not break the statistical rotational symmetry. First we consider the 
case of having a short-ranged repulsive coupling between the two real replicas by which the total Hamiltonian becomes, 

HA+B = Ho[Vo,q^A]+Ho[Vo,q^B]+e [ dz6{Mz) ~ ^b{z)) (19) 



with e > 0. 

This type of perturbation was first considered by Parisi and Virasoro | |47t in the context of spin-glass models in 
order to give a precise definition of spontaneous replica symmetry breaking. It explicitly breaks the RS noted in 
section I. It was also used in the DPRM problem by Parisi in |Q and was further examined by Mezard using the 
numerical transfer matrix method |^ . 

If the two replicas jump into different valleys avoiding to touch with each other, the energy is reduced by an amount 
of order e{L/Lo)- Thus we read off a = 1 and 6U — e so that the overlap length (|l5| ) becomes, 

\ -3/2 



L,(e) -Lo{ — ] (20) 



Again this length scale agrees with the result obtained by Mezard for the same quantity in 1 28 



C. Potential Change, Random 

Now we introduce three other kinds of perturbations 
section ^ this class of perturbations also breaks the RS 
sub-extensive ~ L^/^ (a — 1/2) and much weaker than 
^ L {a = 1). 



Tilt Field and Temperature Change 

which do not break rotational symmetry. As we explain in 
noted in section I. However, the strength of perturbation is 
in the case of explicit repulsive coupling which is extensive 



1. Potential Change 

We consider three different perturbations of a = 1/2. The first one is to introduce a small difference between the 
realizations of the pinning potential for A and B |2^, 30, |3|, |3^. Suppose that A has a certain realization of the 
pinning potential Vq. Then we can construct the potential for B as the sum of Vq and a new statistically independent 
random number Vi. Then the total Hamiltonian becomes, 



Ha+b = Ho[Vo, cPa] + i?o[(K, + SVi)/Vl + S^cbs]- (21) 

Here S is the strength of the perturbation and Vi has the same statistical properties as Vq given in (^ . Namely it has 
zero mean and short-ranged correlations. 



Vi{cl),z)Vi{q^',z) = 2D5i^-(l)')Siz-z'). Vo{(l), z)Vi{^' , z') ^ 0. (22) 



Note that the pinning potential for replica B is normalized by the factor so that it has the same amplitude 

as A replica. 
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The characteristic fluctuation of the extra energy gain along a configuration due to the random variation of the 
potential scales typically as SUq^/L/Lq since it gives contributions with random signs. Thus we read off a = 1/2, 
5U — SUq and the the overlap length ( p^ ) becomes, 

Lc{6) ^ LaS-^ (23) 

This length scale was found by Feigel'man and V. M. Vinokur by essentially the same argument [pO[. Previous 



numerical calculations ||29| , 32 appear consistent with it but the anticipated crossover phenomena had remained to 
be clarified. 



2. Random Tilt Field 



Similarly, we consider the application of a random tilt field to the end-point of B, 

Ha+b=Ho[V„,Q,<I)a]+Ho[Vo,OAb]-5 dzh{z)^^^ (24) 

Here 5 is the strength of the perturbation and h{z) is a Gaussian random number with zero mean and short-ranged 
correlations, 



h{z) = h{z)h{z') = 25{(j) - cf)')5{z ~ z'). (25) 

Within the lattice model we study numerically, the energetic gain of energy typically scales again as SUo^/T/Ui. Thus 
we find a = 1/2 and 5U ~ SUq which gives the overlap length, 

Lci6) ^ LaS-^ (26) 



3. Temperature Change 



All perturbations discussed so far are non-thermal perturbations. Finally we consider the introduction of a slight 
temperature difference between the two real replicas A and B, 

Ta = T + STa 

Tb = T + STb (27) 

where STa ^ STb- Although this perturbation appears to be rather different from the two other cases above, it is 
also expected to give a = 1/2 based on the following observation. 

Fisher and Huse[^ conjectured that valley-to-valley fluctuations of the energy and the entropy are just that of a 
sum of random variables put on a string of length L. Thus the amplitude of valley-to-valley fluctuation scales as, 

AS{L) ^ k^iL/L^f/^ 

AEiL) ^ UoiL/LoY^^ (28) 

However, it is argued that the free-energy is optimized so that these wild fluctuations cancel with each other as much 
as possible in such a way that valley-to-valley fluctuations of the free-energy is much smaller, 

AF{L) = UoiL/Lof with 61 < 1/2. (29) 

In other words, there is a strong negative correlation between the fluctuations of entropy and energy such that 

{AS/kB)iSE/UQ) ^ -(L/Lo) < (30) 

due to the thermodynamic relation AF — AE — AS . Note that a similar argument lies at the heart of the droplet 
theory for spin-glasses which suggests temperature-chaos [|l], |2j . Actually the exponent for the free-energy fluctuations 
is believed to be exactly 9 = 1/3 which is definitely smaller than 1/2. Furthermore, the stronger fluctuation of 
entropy and energy ( [28| ) was confirmed numerically by a transfer-matrix calculation while the smaller fluctuation 
of free-energy with 9 = 1/3 was also observed simultaneously Then under a slight temperature-difference 

between the two replicas A and B, it is possible that one of the replicas jumps into a different valley taking advantage 
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of the large gain in entropy. Such a gain should typically scale as fee \STa - STbUL/LoY/'^ and therefore a = 1/2 and 
SU = k^\5Tji — STb\. From (|l5|), one then finds the overlap length as, 

L.iST)^{!^y (31) 

with ST = STji — STg. This length scale was found by Fisher and Huse in Indeed their transfer matrix calcu- 
lation presented in Q suggests the existence of crossover phenomena. However, details of the scaling properties and 
comparison with the case of the perturbation on potential have remained to be explored. So we try to complete the 
investigation in section 



VI] 



As we summarized above, what is crucial is the role of entropy. In the so called Larkin model | |50{ , in which the 
effect of pinning is modeled by quenched random forces with short ranged correlations, entropy plays very little role 
and free-energy is dominated by energy so that there is no temperature-chaos (sec p^). 

D. Moments of transverse jump distances 

In order to characterize the jump events triggered by the perturbations, it is useful to introduce appropriate 
correlation functions. First, let us introduce the disorder average of the q-th moment of the transverse distance 
between the end points of the two real replicas, 



Bg{L,SU) = {<pA{L)-^BiL))'}. (32) 

It was introduced and studied numerically by Zhang in p9t and continued more in for the case of perturbation 
on the random potential. The following is an extension of the argument by Feigel'man and V. M. Vinokur described 
in m. 

1. Weakly Perturbed Regime 

In the weakly perturbed regime L <C LdSU), a jump event happens with a probability smaller than 1 as given in 
(p^). By a single event, a transverse displacement of order uq^L/Lq)'' will take place. Thus we expect 



Bg{L,SU) 



"o(;^j Piurr.p(.L,SU)^u,iSUy L<.L,{5U), (33) 

where in the last step we have used the scaling relation (||). 

2. Strongly Perturbed Regime 

In the strongly perturbed regime L ^ Lc(SU), the jump events with longitudinal size LdSU) and transverse size 
Uc{SU) will take place with probability one. Let us first consider the behavior of the 1st moment Bi{L,SU) in this 
regime. 

In the strongly perturbed regime, the two replicas A and B are subjected to very different free-energy landscapes. 
In such a situation, we expect that the two replicas A and B will make excursions independently. Thus we expect a 
simple scaling form. 



(</-2)C+a+l 



C 

Lc{SU)J 



BiiL, SU) = UrXSU) [ ^^Tn) ^ » I^ciSU). (34) 



However, the situation is slightly different in the case of uniform tilt field considered in section IV A. Because the 
uniform tilt field continues to increase the separation between A and B systematically as i — s- cx). After making a 
transverse jump of order uo{L / Lq)'' another jump into a further valley in the direction of the field can take place 
if the strength of the field h is increased further. The latter happens when the new increment of the "Zeeman 
energy" 5huo{L/Lo)^ due to another increment of the field 5h becomes again comparable to the typical free-energy 
gap AF*yP(L) given in §), 

Shuo{L/L^f - AF'yP(L). (35) 
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The number of times that such a sequence of jumps happens by increasing the field from to h will be typically h/6h. 
Each jump will have a typical transverse size of order uo{L / Lq)^ . Thus the 1st moment grows as, 

Bi{L,huo) ^ uq ^-^^ ^ ^ Uc{huo) L Lc{huo) (36) 

In the last equation, we used the scaling relation (^). Note that the first moment {q = 1) grows linearly with L not 
only in the strongly perturbed regime but also in the weakly perturbed regime as one can see using a = C in (|3^). 
Actually the linear growth of the 1st moment can be proved rigorously using the statistical rotational (tilt) symmetry 
of the system |2^ , p^ . This is a rather special property of the 1st moment. All other moments are sensitive to the 
crossover from weak to strong perturbation regimes. 

Let us now consider higher moments q > 1. Since jump events are typical in the strongly perturbed regime, we 
generically expect a simple relation between different moments, 

B,iL, SU) = uliSU) (^^^^0^) ' L » L,{5U). (37) 

where note that the natural unit for the q-th moment is now u1{5U). Note that in the weakly perturbed regime such a 
simple relation between different moments does not hold because of the rareness of the jump events. The 1st moment 
obeys a scaling law such that Bi{L,SU)/uc{SU) is a function of L/Lc{5U) also in the strongly perturbed regime as 
we mentioned above. This implies that the higher moments (g > 1) obey a scaling law such that Bq{L^5U)/u%{5U) 
becomes a function of L/Lc{6U){^ 1) in the strongly perturbed regime. 



3. Summary 



To summarize, we expect a generic scaling form for the behavior of the q-th moment including both weakly and 
strongly perturbed regimes as, 

B,{L, 6U) = uli5U)B, [j^^ ■ (38) 

Here the scaling function presents the asymptotic forms in weakly perturbed regime L <^ LdSU) and strongly 
perturbed regime L ^ Lc{SU) which we discussed above. 



E. Overlap function 

Another useful quantity to probe the jump events is the overlap function defined as p8| , 

1 '•^ 



We expect it to scale as. 



q{L,SU) = —l dz6{xo{z) - xsuiz))- (39) 



qmAU)^q[^-^). (40) 



Note that 1 — g is essentially the probability that the string jumps to a different valley. Thus in the weakly perturbed 
regime L <C Lc{SU), we expect that it behaves as, 

l-q{L,5U)^^^^^L,SU)^(^j^y \ (41) 

In the strongly perturbed regime L 3> Lc{6U), we expect that q (i. e. the probability of staying in the same valley) 
decays faster down to as L/Lc{S) oo because the free-energy landscapes of the two replicas are increasingly 
different there. 
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F. Correlation of the free-energy fluctuation 

In order to probe the difference of free-energy landscapes between the perturbed and unperturbed systems, we 
study the correlation function (^) of the sample-to-sample fluctuations of the free-energy between the two systems 
which reads as, 



Cf{L,5U) 



AF{L,0)AF{L,5U) 



(42) 



AF^{L,0)JAF^{L,6U) 



here AF{L, 0) and AF{L, 6U) are deviations from the mean free-energy of the unperturbed and perturbed systems. 
A similar correlation function was studied numerically for the case of perturbation of temperature-shift in |3|. We 
expect it to scale as, 



(43) 



and decay down to as L oo. 

A possible functional form of the correlation function in the weakly perturbed regime L <C Lc{S) can be guessed by 
a simple argument proposed by Bray and Moore for the equivalent problem in a spin-glass model. First, we are 
considering perturbations such that perturbed and unperturbed system have the same statistical properties. Thus we 
must have. 



A^F{L,6U) = JA^F{L,0) - Uq 



L 



(44) 



Suppose that we introduce a perturbation which scales as SU{L/Lo)°' as given in (|l2|). Then the fluctuations of the 
free-energy of the perturbed system have two contributions: the original fluctuation AF{L, 0) plus the change due to 
the perturbation. 



AF(L, 5U) - -i {af{L, 0) + 5U (^|- 



(45) 



Here A/" is a normalization factor which assures that the statistics of the perturbed and unperturbed systems remain 
the same as in(^). It is assumed that the two terms between brackets in expression |45| are uncorrelated. When 
performing the average over the disorder, the cross-terms due to the two terms in (^) cancel out to give the following 
scaling function for the correlation function. 



Cp{L,AU)^^ 



2{a~e)' 



-1/2 



(46) 



For strongly perturbed regime L 3> Lc{S), the correlation hmction may decay faster. 



V. REPLICA BETHE ANSATZ APPROACH 



Now let us take the replica approach introduced in section II to study chaos. We start from the partition function 
of 2 X n replicas: A and B and their n copies. It can be expressed by a path integral over all possible configurations 
of 2 X n replicas labeled by two indices G = A,B and a = 1, . . . ,n, 



/ n 
II []P0G,aexp(-^A+B[0G.a]). 
/-< A D 1 



(47) 



G=A,B Q = l 

where we have introduced the dimension-less effective action. 



B[(PG.a\ = 



dz 



E 

G.a 



2k^T 



d4>G,a{z) 

dz 



D 



G,G',a,p 



(48) 
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To obtain the last equation we have used (^. Here one end of each repHca is fixed as (/>g,q(0) = while the other end 
4>G,a.{L) is allowed to move freely as we noted above. 

The effective action (^) has several important symmetries. First, it has a symmetry under global rotation in the 
(z, 4>) plane. Second, it is symmetric under all possible permutations among the 2 x n replicas. Let us call the latter 
as 'RS' (replica symmetric) for simplicity. As we explained in section II our primary interest is how the RS is broken 
by infinitesimally weak perturbations. 

Now we focus on the study of the disorder averaged partition function ^^^^(L). To this respect we will use the 
well known mapping to an n-body imaginary time quantum mechanical problem in 1-dimensional space, which was 
also firstly noted by Kardar The advantages of this approach is that one can make use of the Bethe 

ansatz which provides us with the exact ground state of the quantum problem. Moreover, from the latter one gets 
many hints about how to construct the relevant excited states. In what follows the main steps in this procedure are 
outlined to emphasize several points which will become relevant in the analysis of the perturbation. The path-integral 
of the partition function defined in (|4^) through the action in ( ^ ) can be reinterpreted as that of a quantum system 
in imaginary time. In the absence of temperature difference between A and B, the Schrodinger equation reads as 



- J^ZI+bH^gA. t) = HoZ^l^siixcc}, t). (49) 
with the following Schrodinger operator for 2 x n-bosons, 

^k^T D ^ 

^o = -L^5^-(fc^ L 5{xo,.-xo.,p) . (50) 

G,Q G,a \ a ! ((G,a),(G',/3)) 

The 1st term represents the kinetic energy. The 2nd term stands for attractive short-ranged interactions between 
the bosons where the sum is taken over all possible pairs of bosons (excluding unphysical self-interactions which are 
absent in lattice models). 

Let us note that here we have two kinds of "bosons" . The bosons of A can be distinguished from those of B and 
vice versa while the bosons cannot be distinguished from each other within the subgroups. However, the Schrodinger 
operator has an even higher bosonic symmetry: it is symmetric under permutations of all the 2 x n replicas. This is 
nothing but the RS we mentioned above. 

By integrating out the coordinates of the free ends of the string ({x^, G}, L) while keeping the other ends fixed at 
(0,0), we formally obtain the disorder- averaged partition fimction of the replicated system as, 



/n 

/n 
n \{dxG.c.<{xG,o}\^t.><i^MQ}> (51) 
/"> /I D ^. 1 



G=A,Ba=l 



where \'ip^ > and are the eigenstates and eigenvalues of the Schrodinger operator Ho defined in (50). In the large 
L (large time) limit, the partition function will be dominated by the eigen states of the Schrodinger operator with 
lowest eigen-values ('energies') including the ground state. 

The ground-state wavefunction is well-known to satisfy the Bethe ansatz reading. 



<^'Rs|{a;G,a}>~exp -A \Xa,G-X(3,G'\\ with X ^ kD / {kBTf , (52) 

V ((G,a),(G',/3)) / 

where the sum is taken over all possible pairs among the 2 x n replicas labeled as (G{— A, B), a{— 1, . . . , n)). The 
index RS stands for the fact that this wave-function has the RS, i. e. permutation symmetry among all 2 x rt-replicas. 
In the following we label this state as replica symmetric (RS). 

In general, the ground state of one-dimensional n-body problems with contact interaction is constructed in the 
following way: the 2 x n particles are ordered and occupy a certain segment within which they are free. The global 
wavefunction consists on the product of 2 x n plane waves whose moments A™ have to fulfill certain matching and 
boundary conditions which in our case result in Am = (2n -t- 1 — 2m) A with m = 1, 2n. The ground-state energy 
is then the sum of the kinetic energy of the 2 x n 'free-particles', 

^. = -^E^™ = -^^'2n(4n^-l) . (53) 

m— 1 
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Although the ground state makes most important contribution to the partition function, it may not be the only 
one. If one only takes into account the contribution of the ground state neglecting all other excited states, one would 
wrongly conclude from ( |53| ) and the relation that only the 1st and 3rd cumulants of the correlation functions 
of free-energy fluctuations exist. This conclusion is definitely unphysical because the 2nd cumulant cannot be zero. 
Such a pathology implies existence of continuum of gap-less excited states which give important contributions to the 
partition function. 

Orland and Bouchaud pointed out that the translational symmetry of the Schrodinger operator allows to 
construct a continuous spectrum of excited states by considering center of mass (CM) motion. Such an excited state 
with wave-vector k has the form, 

< ^RS.k I {xG.a} >= exp(i/c^XG,a) < *Rs|{a;G,a} > (54) 

G,a 

with eigenvalue 

k T 

Ecuik)^ Eg + 2n^k^. (55) 

The resultant partition function obtained by integrating out the continuous spectrum can be put into the following 
scaling form psj , 

lnZj^=-2nf37L + g{2nL'/'') . (56) 

where / in the 1st term represents the average free-energy density. The function g{x) in the 2nd term is analytic 
for small x, implying that the q-th cumulant of the correlation function of free-energy fluctuations scales as L''^^. 
Thus the characteristic exponent for the free-energy fluctuation, which is called stiffness exponent (H), is obtained as 
6 — 1/3, being consistent with extensive numerical results of transfer matrix calculations [2^ and other analytical 
approaches such as the mapping to the noisy Burgers equation [ p7| . 

Parisi p8| pointed out another important spectrum of excited states in which replicas are grouped into clusters of 
bound states. Each cluster is supposed to be described by a Bethe ansatz type wavefunction so that there is replica 
(permutation) symmetry within each cluster. An important assumption is that these clusters are located far enough 
from each other so that their mutual overlap is negligible. The latter is allowed if the transverse size of the system is 
infinitely large. 

In our present context, we have two kinds of bosons corresponding to the two real replicas A and B which can be 
distinguished from each other. Thus it is natural to consider an excited state which consists of separate Bethe type 
clusters < ^'j|g| for A and < 4'^g| for B, with no mutual overlaps, 

<4'rsb| =<*Rs|x <^'fs| <*fsl*RS>=0. (57) 
It's associated energy is readily obtained as, 

EnsB = -'^U'n{n'-l)x2. (58) 

This wavefunction has the reduced replica symmetry mentioned in section |l|, i .e. it is symmetric under permutations 
among A and B groups and the exchange operation A ^ B. We will call this state as replica symmetry broken (RSB) 
state in the following. 

A very important feature is that the gap of the RSB excited state with respect to the RS ground state energy, which 
is of order 0{n^), become vanishingly small in the n ^ limit. Thus such an excited state should be also taken into 
account since we must take n — *■ before L ^ oo in the evaluation of the replicated partition function. Presumably 
each cluster of bound states can have its own center of mass motion. Therefore the RSB excited state should have 
the continuum of excited states of CM motion similar to that associated with the RS ground state mentioned above. 
Then the resultant partition function Z^^^ which will be obtained integrating out these RSB excited states and the 
associated continuum due to CM motions may be put again into the scaling form (^6|). The latter will again yield 
= 1/3. 

To summarize, the replica symmetry is not broken but only in a marginal way. As suggested by Parisi [[48| , the role 
of these RSB excited states will become important if perturbations are considered. In the following we generalize the 
approach of l48| and exploit its implications to study the stability of the frozen phase against various perturbations 



we considered in section IV 
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A. A Perturbative Approach by Replica Scaling Ansatz 



Now we address the situation in which the two real rephcas A and B are under infinitesimally weak perturbations. 
The partition function of the system under such a perturbation can be formally written as, 



yn 



{L)= j n n2?0G,aexp(-5^+B[</'G,a]-<55A+s[0G.a]). (59) 

A D -1 1 



G=A,B Q = l 



where the action SA+B[4'G,a\ is the original one (^ ) which is fully replica symmetric and the 2nd one 5SA+B[4>G.a\ 
is the perturbation term. Suppose that we can map the problem onto the quantum mechanical one such that the 
corresponding Schrodinger operator becomes, 

Ua+b =no + Sn (60) 



where Tio is the original fully replica symmetric 2 x n-boson operator given in (|50|) and STi. corresponds to the SS 
in the path-integral. As we will see in the following, these perturbations try to break the RS present in the original 
system down to the reduced symmetry: replica symmetric only within A and B subgroups. At this stage, the whole 
quantum problem can not be solved exactly. However, we can obtain a useful insight into our problem by perturbation 
analysis proposed by Parisi Esl. 



Here let us note a problem in the case of perturbation by random tilt field considered in section IV C 2 . If one tries 
to obtain a continuous model starting from a lattice model as considered by Kardar , one can find that inter-replica 
coupling terms due to the random tilt field emerge at 2nd order in the transverse hopping rate of the lattice string 
(denoted as 7 in [^). This implies that the mapping in the continuous limit to a Schrodinger equation is invalid in 
this case, because the Schrodinger equation contains only first order time derivatives. Thus we do not consider this 
case in this section. 

From standard perturbation theory we can evaluate the first order corrections to the original ground-state energy 
as, 

where the label A stands for the perturbation strength and (^r.s| is the ground-state wavefunction given in (|5^). The 
1st term corresponds to the ground-state energy given in (|53|). 
Following Parisi, we will consider the RSB excited state (^Tj). 

< ^'rsbKxg.c} >«: exp(-A ^ \xA,a - XA,p\) exp(-A ^ \xB,a - XB.pl) (62) 

a<l3 a<l3 

with A = KD/{k^T)^. This wavefunction has the reduced replica symmetry. At 1st order in perturbation theory, we 
can compute the energy of the RSB excited states as follows, 

^ 3 2k ^ ^ (^-rsbI^-rsb) ^ ^ 

where the 1st term is the energy of the unperturbed system given by (p8|). 

Let us introduce the ratio of the contributions to the partition function Z^^^ due to the RS ground state and the 
RSB excited state, 

D{n, L) = {E^^ - Ef^)L = DQ{n, L) - 5D{n, L) (64) 

where 

Do{n,L) = — \^Tr^L >0 (65) 
is the original 'energy gap' and the correction is due to the 1st order perturbation. 

^ ' V (^-rsbI^-rsb) (^rsI^-rs) J ^ ' 
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In the following we call D{n,L) as "gap". If it is large enough, the contribution of the RSB excited state to the 
partition function becomes negligible. We will find that, in general, the correction term of the gap has the form, 

- 6D{n, L) = -An^L < 0. (67) 

Here the symbol A stands for the strength of the perturbation. Most importantly the correction term —SD{n, L)/L = 
—nPA will turn out to be negative for all the perturbations under consideration. In what follows we will refer to p 
as the order of the perturbation which will play a central role. More precisely the correction to the gap —6D{n, L)/L 
will contain several terms of different powers of n. Here p is the exponent of the term with smallest exponent, which 
becomes most relevant in the n limit. 

If the 1st order correction turns out to give a null contribution, we have to proceed to higher order perturbation 
calculations which is obviously impossible without the complete knowledge of the whole spectrum of excited states. 
Fortunately for all the cases except for the case of the perturbation by uniform tilt field we will find non-zero first 
order corrections. Higher order correction terms will be higher order in A which will be unimportant since we are 
interested in the scaling properties in infinitesimally weak perturbation limit A — > 0. Furthermore it is unlikely that 
the higher order terms are lower orders of n. Thus they will be irrelevant in the n ^ limit. For the case of uniform 
tilt field, we will fortunately find exact RS and RSB bound states of the system which will allow the evaluation of the 
gap D(n, L) also in this situation. 



Now using (g7|) and (|65|) in (|64|) we find, 

D{n,L) ^ Da{n,L) 

with 



n*(A) 



-(3"P)- 



(68) 



A 



l/(3-p) 



From the above result we can generalize the an argument used by Parisi for the explicit repulsive case (p=l) to extract 
the following conclusions. As far as n is integer and the strength of the perturbation A is small, the contribution 
of the RSB state becomes negligible in the thermodynamic limit L oo. However, we have to consider the other 
limiting case: the n ^ limit should be taken before L ^ oo. Now if p < 3, which will turn out to be the case for 
all the perturbations under study, an arbitrarily small perturbation A will induce a level crossing at n*(A) below 
which the contribution of RSB excited state becomes larger than that of the original ground state (RS). The result 
( |68| ) matches perfectly with our definition of chaos (^) since it suggests that the partition function of the total system 
factorizes in the n — s- limit as 



lim ^2+5 = X ^S) ifp<3. (70) 

n/n*— >0 

implying a complete change of the free-energy landscape. 

Now let us further exploit from the above result to find a more physical picture. In the absence of perturbations, 
the logarithm of the replicated partition function has a functional form ( |5^ ) which reads as, lnZ^_|_^ — —l3fL{2n) + 

g{2nL^/^). On the other hand, (68) implies n/n* is another natural variable of the replicated partition function. 



Combining the two, we conjecture the following scaling ansatz, 

In^^;^ + f3]L{2n) - ~g{2nL^'\ n/n*) = ~g{2nL^'\ L/L*) . (71) 
where we introduced a characteristic length L* defined as 

L* - (n*)-3 ^ A-3/(3-p) (72) 

An interesting observation is that n —> limit induces the thermodynamic limit L — > if the variable nL^^^ = x is 
fixed. Then for fixed x we expect, 

g{x,L/L*^0) ~ gi^x) L/L*^l {n/n*^l) 'weak perturbation regime' 

g{x,L/L* ^oo) ~ 2 X g{x) L/L*:^l (n/n* ^ 1) 'strong perturbation regime'. (73) 

The 1st equation means that for small enough lengthscales, the effect of perturbation is small and the partition 
function is essentially the same as that of the unperturbed system of 2 x n-replicas given in (|5^). The 2nd equation 
is the consequence of having to two statistically independent systems in the limit i — > oo. 
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From the above scahng ansatz, it follows that the correlation function of the free-energy fluctuations Cf{L) con- 
sidered in section IVF should have th e scali ng form Cp{L/L*) which goes to as L/L* ^ oo. Similarly the overlap 
function q{L, 5U) considered in section IV E should also have the scaling form q{L/L*) which goes to as L/L* oo. 
Thus the crossover length L* should be identified with the overlap length Lc{5U). The above ansatz implies that the 
de-correlation of the free-energy landscape between A and B takes place as a universal phenomenon whose features 
are classified according to the order of the perturbation p. In the following, we consider the perturbations considered 
in previous real-space scaling argument section ^ specifically one by one based on the replica approach and evaluate 
the correction to the gap ( |67|) explicitly and extract the strength of perturbation A and the order of perturbation p. 
Interestingly enough, we will find that the two approaches give the same overlap length. 

Finally let us comment on how to choose detailed forms of perturbations which we discuss in the following. We 
consider perturbations such that the original symmetry is preserved as much as possible : the 2 x n-replica system 
remains invariant at least under permutation among n-replicas belonging to the same subset A and B and exchange 
A B, i. e. the reduced replica symmetry. 



1. Short-Ranged Repulsive Coupling 

Let us begin with the perturbation which introduces an explicit repulsion term between strings A and B as given 
by the Hamiltonian in (p7|). The corresponding Schrodinger operator for the replicated system can be obviously put 
into the form of ( |60| ) - fully replica symmetric term + perturbation - to obtain, 

Ha+b = Ho + 6n with Sn=-^ y 6{xA,c. - XB,c.) e > 0. (74) 

Q— l,...,n 

Clearly the repulsive perturbing term breaks the original RS Q . 

Computing explicitly the expectation value of the a delta-interaction term with respect to the Bethe ground state 
one obtains Q, 

(*Rs|(5(xA,a - a;B,a)|*RS> _ (2^ + 1) (75) 



(*RS|*RS) 6 

while 

(*RSB|<^(a;A,a - a;B,a)|*RSB 



(*rsb|*rsb) 



0. (76) 



because the bound states of A and B subsets have no overlap < rsb\^ rsb >= 
Thus the correction term to the gap ( |67|) is obtained as. 



SD(n,L) A ^, / e 1 , , 

= --(2n+l)x — — xn. (77) 



L 6 ' ' \kBT 

Note that the reduced replica symmetry: permutation symmetry among n-replicas belonging to the same subset plus 
the exchange symmetry A B is still preserved. Thus the leading order of the perturbation (smallest power of n, 
which becomes most relevant in the n ^ limit) is read off as p = 1 and the strength of perturbation as A '--^ e. 
Finally, using the relation ( |72[) we obtain the crossover length L* ~ e~2 . Remarkably the latter turns out to be the 
same as the overlap length (|20|) found in the real space scaling argument. 



2. Potential Change 

If a slight difference of the random potential is introduced as described in (|2l|), the corresponding Schrodinger 
operator of the 2 x n replica system reads, 

G,a G,a \ a ) ^^^^^ \ a ) ^^^^^ 

^ y ^{xA,a ~ Xb,I3) 



VTTTHkBTY ^^^^^ 

Ho + SH (78) 
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with the symmetry breaking term 



(57i = — - ^ (5 {xA,aiz) - XB,p{z)) . 



(79) 



(",/3) 



Here we are in the infinitesimally weak perturbation Hmit, 5 so higher order terms can be ignored. 

A remarkable feature is that the 2nd term of the last equation, which is the perturbation term STL, is again repulsive. 
Note that the sum is taken over rather than n{n — l)/2 pairs. The expectation value of the delta-function with 
respect to the RS ground state and the RSB excited state has already been computed in ( [75[ ) and (76), hence we 
immediately find the correction to the gap as, 



5D{n,L) 



--(2n+ 1) X — X n"^ 
6^ ^ 2 



(80) 



Note that this perturbation contains the reduced replica symmetry. The latter was made possible by a specific choice 
of the perturbation by introducing the rescaling factor l/Vl + used in (pl|). We can now read off the order of the 
perturbation as p = 2 and the strength of the perturbation as A ~ (5^. Now using (72), we obtain the overlap length 
L* ^ . Indeed, the latter turns out to be the one obtained by the real space scaling argument given in (23). 



3. Temperature Change 

Now two real replicas in the same quenched random potential y((/)(z),z) are subjected to a small temperature 
difference. 

The Schrodinger operator for the 2 x n-replica system with A at temperature T/t and B at temperature T^ is the 
following. 



n 



E 



2k dx\ ^ 



E 



2k dx\ , 



E 



D 



ikBTaXkBTa') 



S{XG,a - XG',f: 



(81) 



((G,a),(G',/3)) 

The RS is apparently lost in the operator. Let us choose the following specific parameters of the perturbation. 



Then we can put the operator in the form, 
with the symmetry breaking terms. 



Ta - 


T- 


hST 




Tb - 


T- 


-5T 




D - 


D 


1 - 3 






n -- 


= ^0 + 


5n, 



sn 



E 

a 

E 



kB6T y 
2k dx\ _ 

kBST 
2k dx'^ , 



2-E 



D 



T ^ ikBTf 



2-E 



D 



T ^ [kBTY 



+2 



ST 
~T 



E 

Q,/3 



D 



5{xA,a - XA,f3) 



S{xB,a - Xb,i 



STY 



^^^^-Six^^^-XBj,) + 0^-j 



(82) 



(83) 



(84) 



In this last equation, we are considering the limit of an infinitesimally weak perturbation ST /T ^ to neglect higher 
order terms. The expectation value of the perturbing operator with respect to the RS ground state is obtained as. 



D 



ST 



T) [kBT]- 



< *Rs|<S(a;A,a - XB,I3)\-^RS > 



a, 13 



>^ D ,^ , ^ STY 2 
--, T7;i2n + 1) X 2 — X n^. 



(85) 
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Here we have used the fact that the ground-state wavefunction is symmetric with respect to the exchange A ^ B 
plus (|75|). Due to the latter, the terms of order 0{ST) cancel out and we are left with the 0{ST^) term. Note also 
that the sum is again taken over n x n pairs of replicas rather than n(n — l)/2. 

On the other hand, the expectation value of the perturbing term with respect the RSB excited state is obtained 
immediately as < ^rsbI'^^I^rsb >= using ( [76| ) and the fact that RSB wavefunction is symmetric with respect to 
the exchange A ^ B and (|5^). Using the above results we find the correction to the gap as, 

6Din,L) X D ,^ ^ f 5tV , 

^ - r(2n+l)x2 — x (86) 



Note that the resultant gap is invariant under the exchange A *^ B which was made possible by the anti-symmetric 
direction of the change of temperature (^2[). From the above results, we read off the order of the perturbation as 
p = 2 and the strength of the perturbation as A ~ (<5r)^. Quite remarkably the latter used in ( |7^ ) again yields the 
crossover length L* ^ {5T)^^ which is the same as the one found by the real space scaling argument (pl|). 



4- Uniform Tilt Field 

Finally we consider to apply a uniform tilt h to one real replica and ~h to the other. The effective action describing 
the uniform field perturbation (|l^) is the following. 



SA+B[(t>G,a] ^ dz 

Jo 

_ h ^ d(f>A,aiz) h d(pB,a{z) 

kBT ^ dz kBT ^ dz 



(87) 



Here not only the full permutation symmetry among the 2 x n replicas but also the global rotational symmetry is lost 
due to the field. Thus the universality of this perturbation should be very different from the ones discussed so far. 
The corresponding Schrodinger operator of the quantum mechanical problem reads as, 

^__v^fcBT_a^ D_ sr^ c( _ . d d 

Note that the first two terms are the original operator Hq given in (^0|) . 

Now let us analyze the change of the RS state (|5|). One can easily see that the 1st order perturbation vanishes 
simply because the total "momentum" of the ground state is zero. On the other hand, one can also easily note that 
when a field is applied, the original wavefunction is no longer an eigenstate. Fortunately, the exact eigenstate can 
be found in this odd situation in which particles belonging to different subsets (A and B) are driven into opposite 
directions. The former Schrodinger operator(|8^) can be rewritten into the fully symmetric form of the original problem 
(^) by shifting the momenta, 

d d h d d h ^^g^ 



Notice that this transformation preserves the commutation relations between conjugated coordinates and moments (i.e. 
[-Q^p- — , XG,a\ = — , XQ^a])- In tcmis of thcsc new coordinates the RS ground state again takes the form of the 

Bethe Ansatz solution of (p2[). And, therefore, the final groundstate can be obtained from Bethe's wave- function by 
undoing the previous shifting of moments, 

* *Rs({a;G,a})exp I > I exp I -— — > XB,a I , (90) 




where 4'Rs({a^G,a}) is the original Bethe ansatz wavefunction for 2 x n replicas given in (|52|). The eigenvalue Eh 
corresponding to this wave function is obtained as 



(91) 
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which does not depend on the ordering of the particles. Here Eq is the original ground-state energy Eg given in (|52^ 
Although the original full permutation symmetry is lost in the wave function (|90|), it still described a sort of bound 
state on 2 X n-particles. So we may refer to it as RS state. In the next subsection, we will discuss the mapping onto 
the Sinai model and the physical meaning will become clearer. The 2nd term of (|9^) gives the change of the eigenvalue 
of the RS state due to the perturbation AEfis = nh"^ / [Kk-QT) . 

Next let us consider the change of the eigenvalue corresponding to the RSB excited state which again is formed by 
two separate bound states for A and B subsets. Here it is useful to note that if all the particles are subjected to the 
common field, the unperturbed single-bound-state wavefunction is still an eigenstate of the operator. Based on this 
observation, one immediately finds that the unperturbed RSB wavefunction is still valid eigenstate under the field 
because of the two fold reasons: i) there is no overlap between A and B and ii) rotational and replica symmetries 
are preserved within the same subsets. Thus the eigenvalue of the RSB state does not change by the perturbation 
AErsb = 0. 

Using the above values of AEfis and AEj^sb we obtain 



D{n,L) 
L 



(92) 



We can now read off p = 1 and A ~ /i^, which yields the overlap length Lc ^ h ^ . Then using ( [72|) we find the same 
overlap length L* ~ being consistent with the result dll) of the real-space scaling argument. 



VI. MAPPING TO A MODIFIED SINAI MODEL 



In the previous section, we found de-correlation the free-energy landscapes of perturbed and unperturbed systems. 
Here we analyze the problem further for the case of uniform tilt field based on the connection between the 1 4- 1 
dimensional DPRM and the statistical mechanics of the Sinai model I 



56 



58 



With this mapping, effective 

1-dimensional energy landscape for the free end x{L) of the 1-1-1 DPRM is obtained as a Sinai potential which is 
generated by a simple random walk in a 1-dimensional space. 

Here we consider this mapping onto the Sinai model in the presence of the uniform tilt field by evaluating the 
partition function (pl|) . First, we evaluate the partition function assuming the RS and using the ground-state wave- 
function given in (BfflT Second, we perform another evaluation assuming replica symmetry breaking (RSB) which only 
allows the reduced replica symmetry and using the 'clustered' wave function (|6^). The former is supposed good for 
the weakly perturbed regime L <C Lc{h) while the latter is good for the strongly perturbed regime L ^ Lc{h). In 
order to interpolate the two limits, we propose a phenomenological model using a bounded Sinai potential. 



5. Replica Symmetric case 

We start by considering the fully replica symmetric (RS) ansatz following Bouchaud and Orland ||5^. The ground- 
state wavefunction under a uniform tilt is given by (|9C|). In order to take into account the motion of the center of 
mass (CM) of the 2 x n replicas, we consider the spectrum of excited states whose wavefunctions are given by, 



^'Rs(/i, k : {xG,a}) ~ ^'Rs({a;G,Q}) exp 



k^T 



exp I ik XG,a 



(93) 



The 1st factor is the Bethe wavefunction given in (^) which describes the unperturbed bound state of 2 x n replicas. 
Now we use the Gaussian transformation introduced by Parisi in Gsl to represent Bethe's wavefunction as follows. 



*Rs({2;G,a}) ^ e ^S((G, = ).(G',3)) 



VV exp 



exp ^V{xG,c,) 



(94) 



,G,a 



The 2nd factor arises from the uniform tilt perturbation: h to subset A and —h to subset B. The last factor is the 
plane wave of wavevector k which accounts for the free CM motion. Here the ground state is included as the fc = 
case. One can easily find the eigenvalues to be. 



E^s{Kk)^Eo 



nh? 
KksT 



knT. 



(95) 
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The 1st term is the original ground-state energy Eg of the unperturbed system given in (|53|), the 2nd term is due to 
the perturbation and the last term is due to the CM motion. 

Let us now suppose that replicas have both ends fixed: one at (0, 0) and the other at xa in the case of subset A and 
xb in the case of subset B. Then the partition function ( |5l| ) is evaluated by integrating out the spectrum of excited 
states as, 



"4 (T ' 

4A V dx 



exp {nV{xA) + nV{xB)) 



X exp((h/kBT){nxA + uxb)) I dkxlLn— — exp [ — Ln— ^ — + ik{nxA + nxB) 



- e-(^°+"'' /(«fcBT))L [g^p ^-nLy^Ens-Sin^iiT, h, V; va, Vb))] ^. (96) 
where [• • -Jy means the average over the effective potential V, 

— ^ = /'pVe-/'^^(^/^')(^)^•• (97) 



and i?R,s-Sinai(r, ft., ua^Vb) is the effective Hamiltonian, 

K (yA + ys)^ ~ , ~ , h h 

i'RS-Sinai^r, /i, V; VA, Vb) = ^^j^ ^ + V(?/a) + V(2/b) - -^Z/A + -j^VB (98) 

in terms of the scaled variables, 

x^L^'^y V = -L^'^V h^L-^'^h. (99) 

By increasing riL^^^, the partition function will be dominated by the minimum of the effective Hamiltonian 
i?RS-Sinai(r, /i, V; j/A, ys)- Then the following physical interpretation can be made: the end point of the strings 
A and B are subjected to the same effective quenched random potential which displays the long-ranged correlations 
in transverse space just as the Sinai model, 



my)-V{y')Y]y^\y-y'\ . (100) 

Furthermore, the CM of the total system is subjected to an effective Hookian spring which tries to bind together 
the two real replicas. The effect of the uniform tilt field amounts to an effective transverse force h applied at the 
endpoints of A and B replicas which tries to drive them into the opposite directions. From (p9|), it can be seen that 
the effective force h increases by increasing the system size L (with fixed h). 



6. Replica Symmetry Broken Case 



In section V A4, we found out that replica symmetry breaking becomes important at L ^ Lc{h) with Lc{h) h 
given in (^). Here we perform the evaluation of the partition function (|^) based on the replica symmetry breaking 
(RSB) ansatz. In this case, we consider a spectrum of excited states whose wavefunctions are given by. 



*RSB(ft, fcA,fcs : {xG,a}) ^ '^R'3.{{xA.,a})'^'Rs{{xB.,a}) G^V 



il^A ^ ^ XA.a 



exp 



ikB ^ a 



XB.o 



(101) 



The first two factors are due to the original wave function of the RSB state ( p^ ) which consists in two clusters of 
bound states. As we noted in VA4, it remains as an eigenstate even under the uniform tilt field since it is assumed 
that these clusters have zero overlap. Moreover, this absence of overlap also allows independent CM motions of A 
and B subsets. The later two factors account for such separate CM motions. The eigenvalues are the following, 



EY{SB[h,kA,B ) = EjiSB +n-(ikAj +n-(ikB) + n— — fc^ + "^-^^ — 



B- 



(102) 
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The 1st term is the original ground-state energy Eyisb of the unperturbed RSB state given in (pSD- The 2nd and 3rd 
terms come from the perturbation. The last two terms are due to the separate CM motions. The partition function 
( ^ ) is evaluated by integrating out the spectrum of excited states as, 



^rsb(0,0|xa,x_b) 



dkA\l Ln ^^^^ 

{Eo + {2n)h^/{2KkBT))L 



iX \ dx J 

dkB\l Ln— — exp 
2ktt 



VVb exp 



d.1- f^V 

A\\ dx J 



exp {nVA{xA) + nVBixs)) 



— 5 — fcn — Ln—likA) + Ln—{ikB) + ikinxA + nxB) 

K 



h 

K 



exp (^-nL^/^EusB-SinmiT, h, Va, Vb^Va, Vb, 



Va,Vb 



(10 



where E^a,B-Sina.i{T , h, Va, Vb; yA, Vb) is the effective Hamiltonian, again in terms of the scaled variables, 



RS-Sinai 



{T,h,VA,VB;yA,yB) 



2kBT 



Va + "^AiyA) - 



kBT 



VA + 



2kBT 



yB + "^B{yB) + -T^yB 

kBi 



(104) 



It is interesting to compare the last result with the replica symmetric (RS) one given in(|9q). Here the two subsets 
A and B are now subjected to independent Hookian springs which try to confine the CM of each subset while the 
total CM was confined in the RS case. Moreover, the two replicas are now subjected to completely independent Sinai 
potentials Va and Vb- The effect of the uniform tilt field again amounts to an effective transverse force h applied at 
the endpoints of A and B replicas which tries to drive them into opposite directions just as in the replica symmetric 



7. Discussion 



In section ^ we conjectured a possible scaling form (|7l| ) of the crossover from the weakly perturbed regime at 
length scales shorter than the overlap length Lc where the RS holds, to the strongly perturbed regime where replica 
symmetry breaking becomes relevant, 



In Z 



A+B 



+ l3fL{2n)^g{2nL^'\LlL,) 



(105) 



Indeed, the partition function based on the RS and RSB ansatz given in (|9^) and ( |103| ) have the expected form; the 
0{n) term which provides the average free-energy (3fL{2n) plus a function which contains the two scaling variables 
nL^/^ and h = L^/^h = (L/ Ldh))^^^ ■ In the last equation we used the relation Ldh) ^ given in (p^). 

Here we have only discussed the two limiting ansatze: RS and RSB. The crossover between the two limits remains 
an open problem. Here let us propose a modified Sinai model which interpolates the limits. We define effective 
Hamiltonian for the endpoints' positions of replicas A and B at a given length reads as follows. 



H = Ha + Hi 



Ha^ 



2kBT 



V{yA) + hyA 



2kBT 



VB + ^iVB) -hyB 



where V{x) is a bounded Sinai potential with correlations. 



{V(x)~V{yW ^C{\x-y\) 



with 



C{u) = y+(l-u)6l(M-l) 



(106) 



(107) 



Here the correlation grows as C{u) — u ioi u < 1 and saturates C{u) = 1 for larger separations m > 1. The latter 
saturation (confined random walk) allows to realize statistically independent Sinai valleys at large separations (RSB) . 
Actually suc h a sat uration of the effective energy landscape was observed numerically in the DPRM by Mezard P8[ . 

In se ction I VII A we analyze the crossover phenomena in detail by a transfer matrix method. Subsequently, in 
section VII B|, we analyze the phenomena numerically using the modified Sinai model defined above and compare the 



result with that obtained in the original DPRM. 



VII. NUMERICAL ANALYSIS 



Now we examine numerically in detail the properties of anomalous response of the DPRM towards various per- 
turbations discussed in the previous sections by transfer matrix calculations. We focus on the anticipated universal 



20 



scaling properties of the crossover from the weakly to strongly perturbed regime across the overlap length which has 
not been clarified in previous numerical studies (see however |28| ). 

Specifically, we consider a lattice model on a two-dimensional lattice of size L xU a,s shown in Fig. |l|. The string of 
length L is directed along the z axis with transverse displacements in the direction of the u-axis. The configuration 
of the string is represented by the positions of the vertices "X" which the configuration goes through, i. e. {u{z),z) 
with z — 1, . . . , L. The "gradient" a{z) — u{z + 1) — u{z) is constrained to take only the values +1 or —1. Note that 
elasticity is realized entropically within this lattice model. The random potential V{u^ z) is defined on each vertex 
(m, z) on which it takes a random value drawn from a uniform distribution between — Vq and Vq. The energy of a 
configuration {u(z)} is given by, 

L 

E[V,u]=^V{u{z),z) (108) 

One end of the configuration is fixed at (0, 0) and the other end is allowed to move freely. On the transverse 
direction we have imposed periodic boundary conditions such that V{u + U,z) = V{u,z). The natural unit for the 
temperature is the scaled thermal energy kBT/Vo where Vq is the unit for the random potential. In the following, the 
Boltzmann's constant is set to fce = 1 and the unit for the random potential to Vb = 1, so that we will often denote 
the scaled thermal energy simply as T. 



z 




FIG. 1: The lattice 1 + 1 dimensional DPRM model. This example has longitudinal size L = 12 and transverse size U = 14. 
The thick zig-zag line is an example of the configuration. The string is directed in the direction of the z-axis with transverse 
displacements in the direction of u-axis. 

First we prepare two real replicas A and B identically except for small perturbations which we will describe in 
detail. Depending on the type of the problem, we use either zero temperature p8| or finite temperature versions ]59|] 
of the transfer matrix method to compute correlation functions. Here and in the following X denotes the average of 
a quantity X over different realizations of the random potential and < X > denotes the thermal average of X (or 
simply the value of X at ground state in the case of zero temperature). We have examined various system sizes up 
to L = 10** and have averaged over Ng = 10^ different realizations of the random potential except for the explicit 
repulsive coupling case for which we used system sizes up to L = 10^ and Ns — 10^. The limitation of the system 
size used for the latter case is that we have to take into account explicitly the inter-real-replica coupling in the 
transfer matrix which requires one to keep track of trajectories of two strings simultaneously and thus much larger 
computational effort psj . 

First we examine the mean-squared transverse displacement of the end point due to the perturbation. 



B2{L) = <UAiL)-UB{L)>^. (109) 

Here ua {L) and u b (L) stand for the position of the end point of A and B replicas respectively. Second we compute 
the exact free-energies (or ground-state energies at zero temperature) of both replicas by the transfer matrix method 
and examine the correlation of the free-energies. 



AF.(X)AF.(L) ^^^^^ 



AFl(L)./AF|(i) 
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where AF is the deviation from the mean free-energy, 

AFa{L) ^ FAiL) -F^iL) AFeiL) = FeiL) -F^. (Ill) 

We also computed the overlap function g(L, SU) defined in ( ^9| ) using the method of However it requires much 
computational effort because one has to keep track of trajectories of two strings simultaneously and computation 
was limited to smaller system sizes L ~ 500. So we do not display the result in the following. We only note that 
anticipated scaling (^) was checked within the limited system sizes. 



A. Uniform Tilt Field 



First we examine the case of the perturbation by a uniform tilt field. For simplify, the temperature is set to zero 
T = 0. The two replicas have exactly the same random potential. The difference is that replica B is subjected to a 
uniform tilt field h which amounts to a force acting just on its end, 

L L 

Ea[V, Ha = 0, ua] = VMz),z) Eb[V, hg = h, ub] = ^ [V{ub{z),z)] - huB^L). (112) 

Z = l 2 = 1 

We have used the T = transfer matrix method and obtained the ground states with various perturbation strengths: 
/i = 0, 0.05, 0.1, 0.2, 0.3, 0.4 for each realization of random potential. 

Le t us begin with the mean-squared transverse displacement of the end point due to perturbation B2 {L) defined in 



2 



(109). B2{L) is expected to grow with increasing size L as 7^1+2/3-5/3 ^j^g weakly perturbed regime (see (^)) and 
in the strongly pertur bed r egime (see ( ^7| ) and (|3^)). Here we used the exponent associated to this perturbation 



as L 

a — C, — 2/?, (see section IV A). The crossover between the two is expected to take place at the overlap length 
Lc ~ given in (|T^). 

In Fig. ^, the data of B2{L) and its scaling plot is shown. For very weak perturbations h = 0.05, the data 
grows almost entirely as i^+^/s except for a short length transient. On the contrary, the data corresponding to a 
strongest perturbation h = 0.4 grows almost entirely as L^, again except for a short length transient. The data for the 
intermediate range of h displays a crossover between the two. Indeed, the scaling plot confirms the expected crossover 



scaling between the two regimes with no adjustable parameters. 




FIG. 2: B2{L) of uniform tilt field case. The data is shown on the left and its scaling plot is shown on the right. Here the 
scaling parameters are Lc{h) = and Uc{h) = Lc(h)''~'^^^ . 

Next let us examine the correlation of the ground-state energies of the perturbed and unperturbed systems through 
( |110| ). In Fig. ||, the data of the correlation function and its scaling plot is shown. The data shows a de-correlation 
of the (free-) energy landscape of the two systems as expected. The scaling plot is obtained again without any 
adjustable parameters. The initial part of the master curve is well fitted by the expected form ( |4^ ) using a = 2/3, 
Cl{F) = 1/(1 + A{L/Lc{h))^'-"-^/^^) with A ~ 2.0. Note that the decay is faster for L/L^{h) > 1. 
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L ULXh) 

FIG. 3: Ce{L) of the uniform tilt field case and its scaling plot with Ldh) = . The fit is Cf(L) — 1/(1 + 
A(L/L,(/i))2(2/3-i/3)) ^-^^ ^ ^ 2.0. 



B. Modified Sinai model 



In section VI we proposed a modified Sinai model as an effective model for the free-ends of DPRM under uniform 



tilt field. Here we study numerically the prop erties of "ground states" of the modified Sinai model and study the mean 
squ ared displacement corresponding to ( 109| ) and the correlation function of the ground state energies corresponding 



to 



110). The effective Hamiltonian given in (106) and (107) at a given length L reads as ]7C 

1 



H = Ha + He 



Ha = 



1 

2l' 



V{xa) 



Hr = 



where V(x) is the modified Sinai potential with correlations, 

(V(a;) - V(?/))2 = u + {u*{L) - u)e{u - u*{L)) 



2L 



with 



V{xb) - h xb 



(113) 



(114) 



First we prepared Sinai potential V{x) on a 1 dimensional lattice u = 1, 2, . . . , i? of size R by generating random 
walks in 1 dimensional space (regarding the 1 dimensional space coordinate as the "time" coordinate for the random 
walk). We generated the bounded Sinai potential by a 1 dimensional random walk confined in a box of size u* . Each 
step of the random walk has variance 1. The same random potential is generated for two replicas A and B. For the 
B replica, we add an extra tilting potential —hu. Then we numerically looked for the "ground states" of A and B 
replicas. We examined various system sizes up to i? = lO'' and used 10^ samples for the disorder averages. 

The 2nd moment of the distance between the minimum is computed for various L and h as. 



B2{L) = K'"(L) - M'^'"(L)]2. (115) 

In Fig. I, the mean-squared displacement is shown together with the scaling plot. In the scaling plot, we included 
the master curve of the equivalent DPRM problem shown in Fig ^. We used the anticipated scaling factors Lc{h) = 
{0.9h)~^ and Rc{h) = 1.2h~'^. The numerical prefactors are chosen such that the master curve of the modified Sinai 
model lies on that of the DPRM problem. 

The correlation function of the fluctuation of ground state energies is computed for various L and h as. 



AEAiL)AEBiL) 



AE\{L) 



AEliL) 



(116) 



where AE{L) is the deviation of a ground-state energy from the mean ground-state energy. In Fig. we show the 
correlation function of the fluctuation of the ground state energy as well as its scaling plot using the anticipated 
scaling variable L/Lc{h). In the plot, we have included the master curve of the equivalent DPRM problem shown in 
Fig I 

It can be seen that the agreement between the modified Sinai model and the original DPRM under uniform tilt 
field is good. We checked that if the original unbounded Sinai potential is used, the agreement becomes very bad for 
large lengthscales. Especially, the correlation function Ce{L) tends to saturate. These results support the picture 
that RSB is needed to account for the de-correlation of energy landscape of DPRM under uniform tilt field. 
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FIG. 4: B2{L) computed by the modified Sinai model and its scaling plot with Rc{h) ~ (0.9/i) ^. In the scaling plot, the 
master curve of DPRM under uniform tilt field plotted vs L/Lc{h) as in Fig. ^ is also included for comparison (black points). 
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FIG. 5; Ce{L) of Sinai model under uniform tilt field and its scaling plot with Rc(h) ~ (0.9/i) ^ . In the scaling plot, the 
master curve of DPRM under uniform tilt field plotted vs L/Lc{h) as in Fig. ^ is also included for comparison (black points). 



C. Explicit Repulsive Coupling 

Next we consider the case of explicit repulsive coupling. The two replicas are at zero temperature, have exactly the 
same random potential and are coupled by an explicit repulsive coupling e, 

L 

E[V, V, e, UA, ub] = ^ {V{uAiz),z) + V{ub{z),z) + e6u^(z),UBiz)) (117) 

Here e > is the strength of the perturbation. Mezard |2^ proposed a transfer matrix method to deal with such a 
coupled system at T > 0. Here we used a T = version of the method and studied the ground states with different 
repulsive couplings e — 0.05,0.07,0.1,0.2,0.3. 

In Fig. H the data of the mean-squared distance of the end points of the two replicas B2{L) is shown together 
with its scaling plot. From the discussion in section IV it is expected to grow with increasing size L as L^' 



in the weakly perturbed regime and L^^^ in the strongly perturbed regime. Here we have used the exponent of the 



perturbation corresponding to the explicit repulsive coupling pertur batio n a = 1 found in IV B (which is related to 



the order of the perturbation p = 2 in the replica analysis in section V A)) . The crossover between both regimes is 
expected to take place at the overlap length Lc ^ e^^/^ given in (|2^). These features are well confirmed by the data 
and the scaling plot. 

In Fig. ^, the correlation of the energies of the two replicas Ea — X]z=i ^("a(-z), z) and Eb — X^z^i ^(""-bI^)? ^) 
is shown together with its scaling plot. The initial part of the master curve matches properly with the expected form 
(M) using a = 1, Cl(F) = 1/(1 + A{L/Lcih))^'^^'^/^^) with A - 0.35. For L/Lc(h) > 1 the decay is much faster. 
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FIG. 6; -02 (i) of the explicit repulsive coupling case and its scaling plot. Here the scaling parameters are Lc{e) = e ^^'^ and 
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FIG. 7: Ce{L) of the explicit repulsive coupling case and its scaling plot. Here the scaling parameter is I/c(e) = e ■ The fit 
is Cf{L) = 1/(1 + A(L/L,(;i))2(i-i/3)) with A ~ 2.0. 



D. Perturbation on temperature, random potential and random tilt field 



Finally we examine the class of perturbations which include temperature-shift, pote ntial change and random tilt 
field. These perturbations are characterized by the exponent a — 1/ 2 found in section IV C| (which is related to the 
order of the perturbation p = 2 in the replica analysis in section V A). Our primary interest here is to clarify whether 
these apparently different perturbations lead indeed to the same universal scaling properties as anticipated by the 
analytical arguments based on the replica-symmetry breaking ansatz. 

1. Potential change 

The Haniiltonian is given as, 



L 

= 5I^("'4(^),^) Ea[V',ub\=Y,V\ub(z),z) 

2 = 1 



(118) 



z=l 



The temperature is set to zero T — Q. First we generate a random potential V{u, z) with random numbers 
drawn from a uniform distribution between —1 and 1. This is the potential for replica A. In order to construct 
the perturbed random potential for replica B, we draw another independent random number U{u, z) from the 
same distribution and define. 



(119) 
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where 6 is the strength of the perturbation. We have used T = transfer matrix method and examined the 
ground states for different strengths of the perturbation 6 = 0.1,0.2,0.3,0.4,0.5,0.6,0.8, 1.0, 1.2. 

2. Random Tilt Field 

The Hamiltonian is given by, 

L 
z=l 

L L-1 

Eb[V, UB,hB] = iVMz),z)) - d J2 hB{z){uB{z + 1) - ub{z))}. (120) 

z=l z=l 

The temperature is set to zero T = 0. The two repHcas have the same random potential V{u, z). The difference 
is that rephca B is subjected to a random tilt field Hb^z) which for each z takes a different random value which is 
drawn from a uniform distribution between —1 and 1. We have used T = transfer matrix method to examine 
the ground states with different random tilt intensities S ~ 0.1,0.2,0.3,0.4,0.5 for each realization of random 
potential. 

3. Temperature-shift 

In the case of temperature perturbation, the Hamiltonian of A and B replicas are exactly the same, 

L L 

Eb[V,ua]=J2^Mz),z) (121) 

z=l z=l 

We have used the finite temperature version of the transfer matrix method. The temperature of replica A is 
set to Ta = 0.1. The temperature of replica B is varied as Tb = Ta + ST with different temperature shifts 
6T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 1.2. 

1. Transverse Jumps 

Let us first examine the mean-squared transverse displacement of the end point B2{L) due to this class of pertur- 
bations. By substituting a = 1/2 in ( ^3| ) we see that B2{L) is expected to grow with increasing size L as l1+(i/2)=3/2 
in the weakly perturbed regime. In the strongly perturbed regime, it should grow as L'^/^ , as discussed in section 
[VP 2| , which is slightly slower than the growth in the weakly perturbed regime. The difference between exponents 
is of only 1/6. The crossover between both regimes is expected to take place at the overlap length Lc ^ LqS^^ as in 
(p3|) , (|3l|) , (^) with S being the strength of the perturbation. 

In Fig. gT'ig. |l^ the data for B2{L) corresponding to the three perturbations are shown together with their scaling 
plots. In the scaling plots we have chosen an adequate numerical prefactor c in Lc{S) = icS)~^ in order that the 
master curves corresponding to the three perturbations lay on the same curve. The resultant master-curves become 
indistinguishable: the expected crossover behavior between weak and strong perturbation regimes is indeed the same 
for the three kinds of apparently different looking perturbations. 

2. De-correlation of fluctuations of free- energies and ground state energies 

In Fig. pl] - p^ , the correlation of the correlation of the ground-state energies or free-energies of the perturbed and 
unperturbed systems are shown together with its scaling plot. In the scaling plots we have used the same numerical 
prefactor c in Lc{S) = (cS)^^ used in the scaling plot of B2{L). As one can see, the master curves for the three 
perturbations merge. The initial part of the master curve fits nicely into the expected form ( |4^ ) using a = 1/2, 
Cl{F) = 1/(1 A(L/Lc(/i))^^^~^/^^) with A - 1.5. One can see that the decay is faster for L/Lc{h) > 1. To sum 
up, the expected crossover behavior from weakly perturbed regime and strongly perturbed regime is indeed the same 
for the three kinds of apparently different perturbations. 
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FIG. 8: B2{L) of the potential perturbation case and its scaling plot with Lc{S) = 5 ^ and Uc{S) = Lc((5) 
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VIII. CONCLUSION 



In this work we have studied the sensitivity of the glassy phase of DPRM against various types of thermal and non- 
thermal perturbations. For the first time, we have obtained very coherent results which strongly support the picture 
anticipated by the phenomenological scaling arguments. As we increase the length scale L at which observations are 
made, there is a crossover from the weakly perturbed regime dominated by rare events (i.e. jumps between neighboring 
free-energy valleys) L ^ Lc{S) to the strongly perturbed regime where these events become typical L ^ Lc{S). This 
means that perturbations become strong at large length scales L/Lc{S) oo such that the configuration can easily 
jump from one valley to another, i.e. it becomes "chaotic" in the sense that the visited landscape is totally different 
from that before. 

In replica space we proposed a new definition of chaos in terms of the global partition function (A-t-B) rather 
than the correlation function itself. There is chaos if in the adequate limits the partition function factorizes, so that 
we have two non-interacting systems. The decorrelation of systems A and B when introducing a perturbation can 
be understood as a concrete example of explicit replica symmetry breaking as proposed by Parisi and Virasoro p7[ . 
Concerning the mapping to the Sinai model, it means that the free-energy landscape of the perturbed DPRM cannot 
be described anymore by a single Sinai potential. Instead, RSB requires the coexistence of statistically independent 
Sinai potentials. 

The key point in our DPRM case is the fact that the RS bound state of the quantum problem is marginally 
stable with respect to RSB as noticed by Parisi [^. Infinitesimally weak perturbations A <C 1 induce small replica 
symmetry breaking terms and induce a symmetry breaking transition from a RS to a RSB state which takes place 
in the n/n*{A) —^ limit for any small but non-zero strength of the perturbation. It turns out that we can read off 
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FIG. 10: B2{L) of the temperature-shift perturbation case and its scaling plot with Lc{5T) = {0A3ST)-^ and Uc{5T) 

lT{st). 



1 

0.8 
0.6 
^ 0.4 
0.2 



o o oo ■ " 

• • 1 



5=0.1 
0.2 
0.3 
0.4 

0.5 
0.6 



10° 



° n □□mm 



□□□CHI] 



10^ 



10" 

L 



10" 



1 

0.8 
0.6 
0.4 
0.2 




10- 





5=0.2 






0.3 






0.4 






0.5 






0.6 


A ^ 




fit - 











10- 



10- 



10"" 10" 

ULXS) 



10" 



10" 



FIG. 11: Ce{L) of the potential perturbation case and its scaling plot with with Lc{5) = {S) 



the overlap length Lc{5) from n*(A). Within the replica space, the perturbations are naturally classified according 
to their order of perturbation p and the symmetries which are left conserved. For each class we have numerically 
verified that, indeed, there are universal scaling functions of correlation functions in terms of L/ Lc{S) describing the 
crossover from the weakly to strongly perturbed regimes. It is notable that the decay of the free-energy fluctuation 
Cf{L) is very slow in all cases we studied. It will not be surprising that one cannot have an impression of "chaos" by 
only making observations within some limited length scales. 

In mean-field models, RSB is always associated with the existence of many pure states ||l5|, which is not the case 
in DPRM in a strict sense. In the DPRM, the mapping to Quantum Mechanics is always possible for an arbitrarily 
large number of dimensions of the transverse space (usually denoted as TV, being A'^ = 1 in our case) in which the 
ground state must contain the bosonic symmetry of the Schrodinger operator p9| . In this sense, RSB has to be weak 
in the DPRM problem p^, i.e. it is a latent feature which only manifests under certain circumstances. In the 
present 1 + 1 model the existence of a 'hidden' RSB excited state with vanishingly small gap with respect to the RS 
ground state in the n — > limit is extremely important. Loosely speaking, the situation is not very far from having 
many pure states. In a really stable RS phase like (like ferromagnetic phases), this phenomena cannot happen. It is 
tempting to speculate that some lessons obtained in the present 1 + 1 dimensional DPRM case may turn out to be 
more general. 

In the present paper, the temperature-chaos is confirmed. Thus the present model serves as a suitable testing 
ground to examine the possible connection between the temperature-chaos and the restart of aging (rejuvenation) 
observed experimentally [|[ ^, |60[ Experimentally, almost complete restart of aging or relaxation take place 

only by slight temperature-changes. Whether this restart of aging can be associated with the temperature-chaos 
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FIG. 13; Cf{L) of the temperature-shift perturbation case and its scaling plot with Lc{5T) — (0.435r) 



remains an interesting open question. Interestingly enough, recent experiments |6l| suggest rejuvenation (chaos) 
effect in random ferromagnetic systems. A candidate to account for the mechanism of the phenomena may be the 
temperature-chaos of pinned domain waUs of the ferromagnets, which is directly related to the present study. 

Another surprise revealed by the temperature-cycling experiments is that initial aging is resumed when the temper- 
ature is cycled back to the initial temperature, giving rise to the so called memory effect. It appears contradictory to 
the temperature chaos effect at first sight. Recently a coarsening model under cycling of target equilibrium states was 
studied [ |62[ . There a hidden dynamical memory by ghost domains was found and a scenario was proposed to explain 
the intriguing coexistence of the rejuvenation and memory effect. In the present context of pinned elastic manifold 
rejuvenation and memory can be easily explained by considering Fourier components of the temporal configuration. 
When the temperature is shifted, Fourier components at wave length larger than the overlap length will be subjected 
to rejuvenation. At time t after the temperature-shift, Fourier components at wave length shorter than L{t) will be 
adopted to the new temperature. Here L{t) is a dynamical length scale over which the system can be equilibrated 
within a given time t. However, Fourier components of even larger wavelength > L(t) remain the same as before the 
temperature-shift. Thus dynamical memory exists at the coarse-grained level of L{t). 

Recalling that relaxational dynamics is extremely slow in glassy systems because of the dominance of the activated 
processes, one has to consider seriously how large time scale is needed to go beyond the overlap length. If it is too 
large, even experimental time scales (typically 10"'^^ — lO""^^ tq where tq ~ 10"^'^ (sec)) may not be sufficient and one 
must look for other mechanisms |63, Sj, ^ to explain the rejuvenation phenomena observed experimentally. Previous 
numerical studies of the relaxational dynamics of the present DPRM model |Q implies that the needed time lies 
within the time window of experiments and numerical simulations for some realistic parameters. More work in this 
direction would certainly be interesting. 
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